#ifndef __Fem_h
#define __Fem_h

#include  "systemvector.h"

/*----------------------------------------------------------*/

class Fem
{
protected:
  double _dx, _chi;

public:
  void initCell(double x1, double x2);
  void initIntegrationPoint(double chi);
  double phi(int i) const;
  double dphi(int i) const;
  void computeFunction(Vector& ulocal, int ic, const SystemVector& uglobal) const;
};

#endif
